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Abstract The Painleve and weak Painleve conjectures have been used widely to identify new in- 
tegrable nonlinear dynamical systems. For a system which passes the Painleve test, the calculation 
of the integrals relies on a variety of methods which are independent from Painleve analysis. The 
present paper proposes an explicit algorithm to build first integrals of a dynamical system, expressed 
as 'quasi-polynomial' functions, from the information provided solely by the Painleve - Laurent series 
solutions of a system of ODEs. Restrictions on the number and form of quasi-monomial terms appear- 
ing in a quasi-polynomial integral are obtained by an application of a theorem by Yoshida (1983). The 
integrals are obtained by a proper balancing of the coefficients in a quasi-polynomial function selected 
as initial ansatz for the integral, so that all dependence on powers of the time r = t — to is eliminated. 
Both right and left Painleve series are useful in the method. Alternatively, the method can be used to 
show the non-existence of a quasi-polynomial first integral. Examples from specific dynamical systems 
are given. 



1 Introduction 



qq ' The present paper deals with autonomous dynamical systems described by ordinary differen- 

tial equations of the form 



Xi = Fi(xi,x 2 ,...,x n ), i = l,2,...,n. (1) 
where the functions F{ are of the form 

N 

F, = ]T OyS°« (2) 
3=1 

with x = (x 1 ,x 2 ...x n ), Qij = (qiji,qij2, ■■■ , qijn), x Qlj = Uk=i and tne exponents q ijk 
are assumed to be rational numbers, i.e., the r.h.s. of Eq.(2) is a sum of quasimonomials 
(Goriely 1992). 

A question of particular interest concerns the existence of first integrals of the system of 
ordinary differential equations (1). A first integral is a function I[xj) satisfying the equality 
dl/dt = V X I ■ x(t) = 0, where x(t) = (xi(t), . . . ,x n (t) is any possible solution of (1). First 
integrals are important because they allow one to constrain the orbits on manifolds of di- 
mensionality lower than n. In particular, a system (1) is called completely integrable if it 
admits n—1 independent and single- valued first integrals. Integrable systems exhibit regular 
dynamics, while the lack of a sufficient number of first integrals very often results in complex, 
chaotic dynamics. 
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The question of the existence of an algorithmic method which can determine all the first 
integrals of (1) is an important open problem in the theory of ordinary differential equations. 
Most relevant to this question are the methods of a) direct search or method of undetermined 
coefficients (e.g. Hietarinta 1983, 1987) b) normal forms and formal integrals (see Arnold 
1985, Haller 1999 and Gorielly (2001) for a review), and c) Lie group symmetry methods (e.g. 
Lakshmanan and Senthil Velan (1992a, b), Marcelli and Nucci (2003)). 

Even more difficult is the question of an algorithmic method probing the integrability of 
Eqs.(l). The most relevant method here is singularity analysis. According to the 'Painleve 
conjecture' (Ablowitz, Ramani and Segur 1980), a system possessing the Painleve property 
should be integrable. The Painleve property means that any global solution of (1) in the 
complex time plane should be free of movable critical points other than poles. According 
to the 'weak-Painleve conjecture (Ramani et al. 1982, Grammaticos et al. 1984, Abenda et 
al. 2001), however, certain types of movable branch points are compatible with integrability. 
Algorithms providing necessary conditions for a system to be Painleve are a) the classical ARS 
test (Ablowitz et al. 1980) b) the perturbative Painleve - Fuchs test (Fordy and Pickering 
1991) which examines the role of negative resonances, and c) the generalized Painleve test 
of Goriely (1992) which introduces coordinate transformations clarifying the nature of the 
singularity structure in the complex time domain. On the other hand, there is no algorithm, 
to the present, determining sufficient conditions for a system to be Painleve. The most 
important obstacle is the search for essential singularities, which are not detectable by any of 
the above Painleve tests. 

The present paper explores the following question: Is it possible to recover the first inte- 
grals of a system by the information provided solely by singularity analysis of its differential 
equations? Our answer is partially affirmative. We cannot circumvent the difficulty con- 
cerning the choice of initial ansatz for the functional form of the integral. The most natural 
choice for quasi-polynomial equations (2) is to consider also a quasi-polynomial ansatz for the 
integral, with undetermined quasi-monomial coefficients. 

This freedom in the initial ansatz notwithstanding, we show in the present paper that 
singularity analysis gives indeed the remaining information needed to recover the integrals. 

First, as proposed by Roekaerts and Schwarz (1987), a theorem by Yoshida (1983) on 
the relation between Kowalevski exponents and weights of weighted-homogeneous integrals, 
can be used to impose restrictions on the degree of the quasi-monomial terms in the integral 
by analysing the resonances found by the Painleve method. Now, Yoshida's theorem for 
weighted-homogeneous integrals is applicable on two conditions: If biT~ Xi is a balance of the 
system (r = t — to is the time around the singularity to), and / is a weighted homogeneous 
integral, the theorem holds if a) VI(bi) is finite, and b) VI(bi) / 0. These conditions shall be 
refered to as 'Yoshida's conditions'. The latter impose a severe restriction in the search for 
first integrals, because one cannot specify in advance whether these conditions are satisfied 
until the integrals are determined. In conclusion, the theorem of Yoshida has only indicative 
power as regards restrictions on the degree of quasimonomial terms in the integral. On the 
other hand, an analysis of the integrable Hamiltonian systems of two degrees of freedom 
with a polynomial potential given by Hiterinta (1983) shows that they all satisfy Yoshida's 
condition when the balance is taken equal to the 'principal' balance, in which all the bi are 
different from zero. Yet, the extent of applicability of this result to other types of systems 
is unknown. In fact, we were able to find also a counterexample concerning a Hamiltonian 
proposed by Holt (1982). 

Assuming a quasi-polynomial functional form of the integral, say <J>(x; ci, C2, cm),&s 
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above, with undetermined quasi-monomial coeficients Q,i = 1,...M, the question now is 
whether it is possible to determine the coefficients q by the series derived via singularity 
analysis. The answer to this question is affirmative. Namely, by the usual Painleve tests, the 
Painleve-type series solutions around movable singularities are first identified 

oo 

X ^ T ) = — Y, b mT m , i = l,...n (3) 

Then, the series (3) is substituted into the quasi-polynomial function <J>(#; c\, ci, cm)- The 
resulting expression is a Puisseux series, i.e., a series in rational powers of r 

oo 

$(x 1 (T),... )% (r);c 1)C2l ..,c M ) = r^ £ d m (b t ; c t )r m ^ (4) 

m=0 

with q,p,m integers. The coefficients d m (bi;ci) are nonlinear functions of the coefficients bi 
(determined by (3), i.e., by singularity analysis), and linear functions of the undetermined 
coefficients Cj. But the function $ is an integral of the system if it is constant along all the 
solutions of the system, including (3). As a result, the functions d rn satisfy the set of linear 
equations 

d m (bf, Ci) = 0, i = 0,1,..., i^-q (5) 

This is an infinite number of homogeneous linear equations with a finite number of unknowns 
(the coefficients Cj). If M is the number of unknown coefficients, Eq.(5) can be written as 

A(k) ■ C = (6) 

where C = (c\, . . . ,cm) T and A(bi) is a matrix with M columns and an infinite number of 
lines. The entries of A depend only on the coefficients bi which were previously determined 
by singularity analysis. In this representation, the first integrals are functions with quasi- 
monomial coefficients given by the basis vectors of ker(A) (or linear combinations of them). 
In computer algebraic implementations of the method, we work on a sub-matrix Af defined 
by a finite number of lines, which is equal or larger than M. A basis for the subspace ker(Af) 
is determined by the singular value decomposition algorithm. Then, it is checked with direct 
differentiation that the resulting expression is an integral. This completes the determination 
of all quasi-polynomial first integrals for the given system. 

This method was implemented in a number of examples presented below. Following some 
preliminary notions exposed in section 2, the results are presented in section 3, along with 
various details and implications in the implementation of the algorithm. Section 4 summarizes 
the main conclusions of the present study. 



2 Preliminary notions 

Following Yoshida (1983), a system of the form (1) is called scale-invariant if the equations 
remain invariant under the scale transformation Xi — ► a x% %i, t — > a~ x t for some \. Then, the 
system (1) has exact special solutions of the form 
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where bj, i = 1, . . . 



n is any of the sets of roots of the system of algebraic equations 



bn) + = 



and the time r = t — to is considered around any movable singularity to i n the complex time 
plane. Any solution of the form (7) is called a 'balance'. If the system (1) is scale-invariant 
and the functions Fj are of the form (2), the exponents Aj associated with any of the balances 
are rational numbers. It should be stressed that the definition above does not require that 
all the b{ be non-zero. Balances of the form Xi ~ 0/r Ai , for some i, are also considered. The 
latter remark is essential in order to avoid the confusion which is sometimes made between 
'balances' and 'dominant terms' in the ARS test (see e.g. the discussion between Steeb et 
al. (1987) and Ramani et al. (1988)), and the associated difference between 'Kowalevski 
exponents' and 'resonances'. In the standard ARS algorithm (Ablowitz et al. 1980) the 
solutions (7) arise by the definition of the dominant behaviors, i.e., which is the first step in 
the implementation of the algorithm. 

The second step in the ARS algorithm is to look for series solutions that we call 'Painleve 
series'. These are expansions of the form (3) starting with dominant terms of the form (7). 
They are Laurent (Taylor) series when the Aj's are integers (positive integers), otherwise 
they can be series in rational powers of r, which are called 'Puisseux series'. To build up 
the series, one first specifies the resonances, i.e. the values of r for which the coefficients of 
the terms r r ~ Xi in the series are arbitrary. In the case when the coefficients of the balance 
bi are all non-zero, the resonances are equal to the eigenvalues of the Kowalevski matrix 
Kij = (dFi/dxj + Sij\i)\ Xi= b x which are called 'Kowalevski exponents'. If, however, some of 
the 6jS are equal to zero, then the resonances are not equal one to one to the Kowalevski 
exponents, but some resonances differ from the corresponding Kowalevski exponents by a 
quantity equal to the difference between the exponent Aj in the balance and the exponent of 
the first non-zero dominant term in the Laurent-Puisseux series of Xj(r), as specified in the 
first step of the ARS algorithm (Ramani et al. 1988). 

At this point, we are not interested in whether the system passes the generalized Painleve 
(or weak Painleve) test. This means that we do not require that all the solutions of the system 
(1) can be written locally (around a movable singularity) in the form (3), or that there is at 
least one solution of the form (3) which contains n arbitrary constants (including to). On the 
other hand, we do check the compatibility conditions to ensure that no logarithms enter in 
the series. As regards positive resonances, compatibility is fulfilled automatically for scale- 
invariant systems. In summary, we are interested only that the system have special solutions 
of the form (3), but no other claim on it being Painleve or not is required. Thus, the results 
are valid also for partially integrable systems, i.e., systems with a number of first integrals 
smaller than n. 

Let us now assume that (1) possesses a weighted - homogeneous first integral <E> of weight 
M, i.e. an integral function $ which satisfies the relation: 



for some M. Then we have the following 

Theorem 1 (Yoshida 1983): //, for a particular balance (7) the following conditions hold: 
a)VQ(bi) is finite, and b) V&(bi) / 0, then M is equal to one of the Kowalevski exponents 
associated with that balance. 




xi,a 2 x 2 ,-.,a n x n ) =a $(xi, x 2 , x n ) 



(8) 
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There has been a number of theorems in the literature linking Kowalevski exponents with 
the weights of homogeneous or quasi-homogeneous integrals of a system. Some characteris- 
tic papers on this subject are Llibre and Zhang (2002), Tsygvinsev (2001), Goriely (1996) 
and Furta (1996). However, the application of Yoshida's original theorem in the search for 
weighted-homogeneous integrals seems to be the most practical. Furthermore, Yoshida's the- 
orem can be reformulated in an interesting way: consider the Painlene series starting with 
one of the balances (7). It follows that the series can be written as: 



x i = x iE + Ri = \ + - + A i T r ~^+R i (9) 



where r is the maximum Kowalevski exponent associated with this balance and Ai is the 
corresponding arbitrary coefficient entering in the Painleve series (9) for the variable X{. The 
sum XiE = -r~ + ••• + AiT r ~ Xi will be called essential part of the series and the remaining 
part Ri remainder of the series. The remainder Ri starts with terms of degree r — Aj + 1/p 
where p is the denominator of Aj written as a rational Aj = q/p with q,p coprime integers. A 
quasi-polynomial integral $ has the form 



$ = E C ^2,...,teIl^ (io) 



i=l 



where the exponents ki are rational numbers. If the integral $ is weighted-homogeneous of 
weight M, we have 



E 

i=i 



h\ = M (ii) 



due to (8). Taking into account the fact that the remainder Ri starts with terms of degree 
0(r r ~ Ai+1 / p ), it follows that the contribution of Ri in x^ is in terms of degree <9( r - A » fc »+ r + 1 /p) 
or higher. This means that if the series (9) is substituted in the integral &(xi), the contribu- 
tion of Ri in §( Xi ) is of degree 0( r - 2~> <fc<+r+1 /P) = ( r -M+r+i/p) or higher. But r > M. 
Thus, the remainder Ri contributes only to terms of positive degree in r. On the other hand, 
since <E> is an integral, the time r as a denominator must be eliminated in <3?(xj). But since 
there are no negative powers of r generated in $ by Ri, it follows that all the negative powers 
of r are already eliminated by substituting the expression XiE alone into <£. Hence, we have 
the following 

Proposition 2: // the conditions of Yoshida's theorem hold for a weighted-homogeneous 
integral of (1) and a particular balance (7), then the expression <S>{xie), where XiE are the 
essential parts of Painleve series Xi(r) initiated with the same balance, does not contain sin- 
gular terms in r. 

Consider next the case when the functions Fi in (1) are not homogeneous. By the restrictions 
imposed by (2), it follows that the functions Fi can be decomposed in sums of the form 

p(m i0 ) _|_ p(m i0 +l/p) _|_ _|_ p(m i0 +q/p) 



p. _ py'Ho) _|_ py>Ho-ti/P) _|_ _|_ py'Ho^-q/P) ^2) 

where the functions F^ are homogeneous of degree j, with j, mio rational, p, q integer, and 
p is the denominator in the simplest fraction giving TOjo- In this case, if the system (1) has 
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Painleve type solutions, the dominant behaviors and associated resonances of these solutions 
are determined by the homogeneous term of the highest degree 

hand, the functions j < + q/p must have a special form to ensure that compatibility 
conditions are fulfilled and the series solution is of the Painleve type. Finally, as regards 
potential first integrals, the assumption that they consist of a sum of quasi-monomial terms 
implies that they can also be written as sums of the form (12). The selection of terms in 
the quasi-polynomial integral can be determined by Yoshida's theorem (or proposition 2) 
implemented in the scale-invariant systems 

x i = F^ mi0 \x i ) (13) 

and 

±i = F (m i0+q / P ) {xi) (14) 

respectively (Nakagawa 2002). 



3 Explicit construction of integrals with quasi-monomial terms 

3.1 An elementary example 

Consider the two-dimensional nonlinear system 

X\ = X2, X2 = —x\ — 3x\ (15) 
The only first integral of this system 

$ = x\ + x\ + 2x\ (16) 

can be recovered by elementary means. However, we will use this example to illustrate the 
steps used by the present method. The corresponding homogeneous system containing the 
terms of maximum degree 

x\ = x 2 , X2 = -3xf (17) 

has the unique balance x\ = — 2/r 2 , X2 = 4/r 3 , with Kowalevski exponents (equal to reso- 
nances) — 1 and 6. Combatibility conditions are fulfilled for the system (15), which admits 
the Laurent series solution 

- ^-\-^ + ^ > + 0{ ^ 

x 2 {t) = A-^r + 4a 4 T 3 + 0(T 5 ) (18) 

T° DO 

where 04 is an arbitrary parameter. 

Following Yoshida's theorem, we shall look for an integral of the system (15) by requesting 
that this integral be a sum of weighted- homogeneous functions of weight not higher than 
M = 6. Since the equations (15) are polynomial, the integral will also be assumed polynomial. 
According to the definition of the weighted-homogeneous functions (8), the undetermined 
integral contains terms of the form x\ 1 xf 2 the exponents of which are restricted by the relation 
2gi + 3<72 < 6. This leaves only six possibilities, namely (q\ = l,q2 = 0), (q± =0,^2 = 1), 
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(qi = 2,q 2 = 0), (q 1 = l,q 2 = 1), (qi = 0, q 2 = 2), (q 1 = 3,q 2 = 0). Thus the integral is 
assumed to have the form 



$ = C W Xl + C lX 2 + c 20 xj + C U XlX 2 + c 02 x 2 + c 30 xf 



(19) 



Up to now, the steps are exactly as proposed by Roekaerts and Schwarz (1987). At this point, 
however, we do not proceed by the 'direct method'; instead, the series (18) is substituted into 
(19). Then, terms of equal power in r are separated and their respective coefficients are set 
equal to zero. We must determine at least as many equations as the number of unknown 
coefficients Cy, i.e., six equations. These are: 

Order 0(l/r 6 ) : 



16c 



Order 0(l/r 5 ) : 
Order 0(l/r 4 ) : 

Order 0(l/r 3 ) : 

Order 0(l/r 2 ) : 

Order 0(l/r) : 
These equations can be written in matrix form 



8c 30 = 
- 8cn = 
4c 2 o - 2c 30 = 



4c i - -en 







2cio + -c 2 o 



IT 02 



15 C3 ° 





= 

















16 


-8 




/ Cio ^ 




( \ 













-8 










coi 















4 








-2 




C20 












4 





-2/3 










Cll 









-2 





2/3 





-2/15 


-4/15 




C02 







V 




















/ 


V c 30 J 




V o J 



(20) 



or simply 



A r C = 



(21) 



where C is a six-dimensional vector and Af is a 6 x 6 matrix with constant entries. 
The singular value decomposition of Af yields a one-dimensional null space: 

ker(A f ) = A(0,0, 1,0, 1,2,0) (22) 

The basis vector (0, 0, 1, 0, 1, 2) corresponds to the first integral $ = x\ + x\ + 2x\ . 
A few remarks are here in order: 

a) the last line of Af has only zero entries, since it corresponds to the identity = for 
the 0(1/ t) terms. This is not a problem, because lines with zero elements are allowed by the 
singular value decomposition algorithm which determines the subspace ker(Af). 

b) The entries of Af are constant numbers which depend only on the coefficients of the 
Laurent series (18). This is the crucial remark; it implies that the information on the first 
integral is contained in the Painleve series built by singularity analysis. 

c) The arbitrary parameter 04 in (18) does not appear in Af. This phenomenon is not 
generic. In general, all the arbitrary parameters of the Painleve series appear in Af. In the 
computer implementation of the algorithm, we proceed by giving fixed values to the arbitrary 
parameters. Although the choice of values affects the convergence of the Painleve series, it 
does not influence the present algorithm which is based only in the formal properties of the 
series. 
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3.2 Further examples 



Of particular interest in nonlinear dynamics are autonomous Hamiltonian systems of two 
degrees of freedom of the form 



where V(x, y) is of the form (2). The easiest examples are systems with a polynomial potential 
(e.g. Hietarinta 1983, 1987). For example: 



This system passes the Painleve test and it is integrable (Bountis et al., 1982). Keeping 
only the highest order terms of the potential {—x 2 y — 2y 3 ) yields the principal balance 
x = 6i/r 2 ,y = 3/r 2 with resonances r = —3,-1,6,8, which are equal to the correspond- 
ing Kowalevski exponents. The Painleve series generated by the Hamiltonian (24) and the 
above principal balance satisfies the compatibility conditions of the ARS test. Assuming now 
a polynomial first integral <E> of (24), only the monomial terms of weight less or equal to 8 
will be included to it. Since the leading terms of the momenta are p x ~ p y ~ 0(l/r 3 ), the 
selected monomial terms are: 



It should be stressed that there are several other restrictions that reduce the number of eligible 
terms. For example, the integral <I> is either even or odd in the momenta (Nakagawa and 
Yoshida 2001). Furthermore, the linear terms can be omitted by appropriate transformations. 
However, in the practical implementation of the algorithm these restrictions only introduce 
a complication, because, if an integral $ exists, the singular value decomposition algorithm 
selects the basis for the corresponding null space of the matrix Af without needing any extra 
information on restrictions which are specific for the system under study, e.g. hamiltonian or 
other. 

Following the selection of monomial terms, the algorithm proceeds in building the ho- 
mogeneous system (6) as well as the finite restriction Af of the matrix A. In this case, the 
dimension of Af should be M x 38, with M > 38, since there are 38 unknown coefficients 
of the above monomial terms. At this point, it does not matter which balance and Laurent- 
generated series are used to build the matrix Af. In the above example, the principal balance 
leads to a 'special solution', since there are only two positive resonances (r = 6,8) mean- 
ing that there are three arbitrary parameters in total entering in the series. On the other 
hand, the general Laurent series solution with four arbitrary constants is given by a different 



H=\{pl+P 2 y ) + V{x,y) 



(23) 




(24) 



x 4 , x 3 y, x 2 y 2 , xy 3 , y A , p 2 x x, p x p y x, p 2 y x, p 2 x y, p x p y y, p 2 y (weight 8) 
p x x 2 , p x xy, p x y 2 , p y x 2 , p y xy, p y y 2 , (weight 7) 

3 2 232 232 2 3 / • r>\ 

x , x y, xy , y , p x , p x p y , p y ,x , x y, xy ,y (weight 6) 
P x x, p x y, p y x, p y y, (weight 5) 
x 2 , xy, y 2 , p x , p y , x, y (weights 4,3,2) 
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balance, namely x = 0/r 2 ,y = 1/r 2 so that x starts with dominant terms of 0(l/r), i.e. 
x(r) = — + air + Br 2 + a 3 r 3 + a 4 r 4 a 5 T 5 + (9(r 6 ) 

y(r) = + b + 6 2 r 2 + 5 3 r 3 + Cr 4 + 0(r 5 ) (25) 
P*(t) = ^ + ai + 25r + 3a 3 r 2 + 4a 4 r 3 5a 5 r 4 + 0(r 5 ) 



r 2 
2 



y(r) = -^ + 26 2 r + 36 3 r 2 + 4Cr 3 + 0(r 4 ) 



T 

where .A, C together with to = t — T are arbitrary, and 6 = (1- A 2 )/12, ai = A(l-2&o)/2, 
fr 2 = (6 -66g-2Aai)/10, a 3 = (26 ai-ai+2A6 2 )/4, 6 3 = -AB/3, a 4 = (2A6 3 - J B-2 J B& )/10, 
as = (2AC + 2ai& 2 + 260^3 — a 3 )/18. In this solution, the resonances r = —1,0,3,6 do not 
coincide one by one to the Kowalevski exponents (—1,1,4,6) deduced by the Kowalevski 
matrix associated with the above balance. Nevertheless, the general solution (25), as well as 
any other solution work equally well in determining the matrix Af. In our case, by performing 
the singular value decomposition of Af, the subspace ker(Af) yields two independent integrals 
in involution, with coefficients (up to the computer precision) 

$1 = 0.211944988455(y 2 + p 2 y - 4y 3 ) + 0.049612730813(x 2 + p 2 x ) 

+0.216443010189(3%^ - p 2 x y - x 2 y 2 - ^-x 4 ) (26) 

-0.207446966722x 2 y 

$ 2 = 0.035960121414(y 2 + p 2 y - 4y 3 ) + 0.342416478352(x 2 + p 2 x ) 

-0.408608475918(xp a; pj / - p 2 x y - x 2 y 2 - ^x 4 ) (27) 
-0.480528718746x 2 y 
in terms of which we can express the Hamiltonian 

H = 2.1645645936295^1 + 1.146586289822$ 2 (28) 
and find a second integral orthogonal to the hamiltonian 

h = -2.1645645936295$ 2 + 1.146586289822*1 (29) 

The integral J 2 is a linear combination of the hamiltonian and of the integral lb given by 
Bountis et al. (1982) 

h = 0.1651752123636227(2# - y/ 6 ) (30) 

Let us note that an alternative way to obtain the matrix Af is by considering only the 
essential parts of the Laurent series (25), for as many different sets of values of the arbitrary 
parameters as requested in order to have a complete determination of the M x 38 elements of 
Af, with M > 38. This approach is preferable in computer implementations of the algorithm, 
because one does not need to calculate the terms of the Laurent series (25) beyond the highest 
positive resonance. Furthermore, Proposition 2, instead of Yoshida's theorem, can be used to 
select the quasiminomial basis set. 
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A final remark concerns the applicability of Yoshida's condition V$(6j) / 0. We checked 
whether this condition is satisfied in Hientarinta's table (1983) of integrable Hamiltonian 
systems of two degrees of freedom with a polynomial homogeneous potential. There are six 
non-trivial cases: 

(a) V(x, y) = 2x 3 + xy 2 
with integral <E> = yp x p y — xp 2 + x 2 X2 + \y 4 

(b) V(x,y) = —x 3 + xy 2 

with integral $ = p 4 . + Axy 2 p 2 y - \y 3 p x p y - |x 2 y 4 - |y 6 

(c) V(x,y) = 2x 3 + xy 2 + i— y 3 

with integral $ = p* + ^jfPxP% + -^W 3 pI ~ ( 2 V 3 + 2iyJ?>xy 2 )p x p y + (4i^/3x 2 y + 2i v / 3y 3 + 
kxy 2 )p 2 y + -^ix 3 y 3 + ^xy 5 - x 2 y A - |y 6 

(d) V (x,y) = ^ + x 2 y 2 + ^ 
with integral $ = yp x p y — xp 2 + {\2x 3 + xy 2 )y 2 

(e) V(x,y) = ^x A + x 2 y 2 + ±y 4 

with integral <$>=p A y + \y 4 p 2 x - \xy 3 p x p y + (Ax 2 y 2 + fy 4 )^ 2 , + |(y 8 + 4x 2 y6 + 4x 4 y 4 ) and 

(f) V(x, y) = x 5 + x 3 y 2 + j^xy 4 

with integral = yp x p y - xp 2 + \x 4 y 2 + \x 2 y 4 + -^y 6 . 

In all these cases, the principal balance, with all bi different from zero, satisfies Yoshida's 
condition. Note that case (f) the principal balance leads to 'weak-Painleve' solutions, but 
the associated integral $ is easily recoverable by the present algorithm. In fact, the role of 
Yoshida's condition is to exclude from Theorem 1 integrals which are composite functions of a 
simpler integral. For example, consider the integral / = <£ 2 , where <3? is the polynomial integral 
in any of the above six Hamiltonian systems. In all of them is weighted-homogeneous for 
some weight M. Thus / is weighted-homogeneous of weight 2M. Substituting the special 
solution (7) in the integral <!> yields <&{bi/r Xi ) = T~ M <fr(bi). Since <!> is an integral, it should 
be time-independent, thus &(l>i) = 0. Similarily, I(bi) = 0. However, while V<3?(6i) ^ 0, we 
have VI = 2<E>V < 1?, thus VI{bi) = 0. Thus, while <I> satisfies the Yoshida's condition, / does 
not. This ensures that while M is necessarily a Kowalevski exponent, the multiples of it, e.g. 
2M are not necessarily Kowalevski exponents. Viewed under this context, it appears that the 
condition V<£(6j) 7^ makes a 'natural' choice of the simplest integral among an infinity of 
possible integrals which are composite functions of <£. 

However, there are interesting counterexamples which challenge this point of view. One 
example is the homogeneous limit of the Holt (1982) Hamiltonian 

H= l -(pl+pl)-\x 4 l 3 -y 2 x- 2 l 3 (31) 
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In this case, the potential is a homogeneous sum of quasi-monomials, of degree 4/3. The 
second integral reads: 

$ = 2y 3 + 3yx 2 - 6yy 2 X - 2 / 3 + 9yx 4 / 3 - I8±yx 1/3 (32) 
and it is weighted-homogeneous of degree M = —6. The principal balance is 

with resonances (=Kowalevski exponents) r = —1,-2,— 3, —4. This case is remarkable be- 
cause a) all the resonances are negative, and b) they are not equal to the weight of the integral 
M = — 6. Thus, Yoshida's theorem does not apply in the case of this balance. Substituting 
the balance b = ((1/3) 3 / 2 , l/3>/2, 3(l/3) 3 / 2 , l/\/2), i.e., Eq.(33) in the gradient of the integral 
(32), 

V$(x, y, x, y) = (Ayy 2 x~ 5/3 + 12yx 1/3 - Qxyx~ 2/3 , (34) 
-l2yyx- 2/3 - I8xx 1/3 , 6±y - I8yx 1/3 , 6y 2 + 3x 2 - 6y 2 x~ 2/3 + 9x 4/3 ) 

yields V<£(&) = 0. Thus none of the resonances has to be equal to the weight of the integral. 
However, the form of <E> or V<3? does not suggest that these functions are composite functions 
of some simpler integral. On the other hand, there is a second balance of the Hamiltonian 
(31), namely b = ((1/6) 3 / 2 , 0, 3(l/6) 3 / 2 , 0), which corresponds to the dominant behavior 

/i \ 3/2 

*=(gj t 3 , y = Ar 4 (35) 

with A arbitrary. The resonances here are r = 0,-1,-4 and —7, but the Kowalevski ex- 
ponents are rx = 1, —1, —4, —6 (two of them differ by one from the respective resonances). 
Now, the Kowalevski exponent = — 6 is equal to the weight of the integral. If we look at 
V$, Eq.(34), we see that the components d<&/dy and d^/dy contain terms independent of y 
and y. Thus, for this particular balance V$(b) / 0, i.e., Yoshida's condition is satisfied. 
Thus, in the absence of a counterexample, we formulate the following 

Conjecture 3: In any scale-invariant system of the form (1,2), which possesses a weighted- 
homogeneous first integral at least one of the balances b satisfies the condition V<3?(6) / 0. 

Returning to the Holt Hamiltonian, the next step is the selection of quasi-monomial terms 
in the initial ansatz for a quasi-polynomial integral. Guided by the form of the Hamiltonian, 
natural exponents are adopted for the powers to which the momenta p x , p y , and of the vari- 
able y are raised, while the variable x is considered as raised to powers m/3, where m is 
integer (positive or negative). Even under these restrictions, there is an infinity of possible 
quasi-monomial terms of weight -6. For example, the terms pix m / 3 y( 2 ~ m ^ 3 , where pi is either 
p x or p y , are of weight —6 for all m G Z. Thus an arbitrary lower limit has to be set to m. 
This is chosen as the lowest bound of m in the Hamiltonian, namely m = —2. Nevertheless, 
failure to find an integral with these restrictions on the quasi-monomial terms does not im- 
ply that an integral does not exist, because of the arbitrariness with respect to the lowest 
bound of negative exponents considered. This problem does not exist when there are only 



11 



positive exponents present in the quasi-monomial terms of the equations of motion (or of the 
Hamiltonian) . 

With the above restrictions, the quasi-monomial terms considered are: 

Pi, P 2 x Py, PxP 2 y , pj, p 2 x x 2/3 , plx~ 1/s y 

2/3 — 1/3 2 2/3 2 —1/3 

PxPyX 1 , p x pyx 1 y, p y x ' , p y x 1 y 
p x x 4/3 , p x x l/3 y, p x x~ 2/3 y 2 , p y x i/3 
P y x 1/3 y, p y x~ 2/3 y 2 , x 2 ,x y ,y 2 

The final step is to build the matrix Af as in the previous examples, i.e. by replacing the 
Painleve series in the initial ansatz for the integral. An interesting point is that in the case of 
the Holt Hamiltonian (31) we have to consider left Painleve series (Pickering 1996), i.e., series 
in descending powers of r. This is because the balances ~ r 3 do not imply singular behavior 
as t — > 0. In this case, the limit |r| — ► oo represents a singularity, but the left Painleve series 
are convergent for all r with |r| > e for some real positive e. The series are constructed as: 

oo 

x t (r)=r x J2hr- k (36) 

fc=0 

where A > 0. Resonances and compatibilities are checked in the same way as in the usual 
Painleve test. By this method, we were able to obtain the integral (32) by a proper balancing 
of the quasi-monomial coefficients Cj so as to eliminate the coefficients of the terms of successive 
descending powers of r in the integral expression. 

As a final remark, it should be stressed that the selection of a quasi-polynomial ansatz for 
the integrals of a system of the form (1), with the functions (2) being quasi-polynomial, is not 
exhaustive. This can be easily exemplified in a case with polynomial functions. The Bogoy- 
avlensky - Volterra B-type systems are given in normalized coordinates U{ by the following 
set of autonomous nonlinear ODEs: 

■di = UiUi+i - UiUi-i, i = 2, ...,n-l (37) 

U n — U n U n —\ 

the r.h.s. of Eqs.(37) are homogeneous functions of second degree in the variables U{. The 
system (37)) admits balances of the form of the form: 

fH = ^ (38) 
r 

where t = t — to is the time near a singularity to i n the complex t-plane. In the case of the 
principal balance, the a\ are non-zero solutions of the set of algebraic equations: 

-a\ = a\ + aia 2 

-di = Oidi+i — didi-i, i = 2, ...,n — 1 (39) 

Q>n = Q-n&n—l 
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and they are given by the recursion formulas 



ak+2 = a k -l, ai = (-l) n [ 



n + 1 , 



a 2 



— 1 — a\ 



(40) 



for k = 1, n. 

The resonances of the principal balance (=Kowalevski exponents) are given by the char- 
acteristic equation, i.e., setting the determinant of the Kowalevski matrix equal to zero. The 
determinant has a tridiagonal form, i.e., 



del 



V 



ai 

a2 — r a2 
— 03 — r 03 















\ 





-a n _i -r a n _i 
-a n -r ) 



= 



(41) 



which can be solved easily yielding the resonances 

r k = (-l) k k, k = l,...,n 



(42) 



Assuming that the conditions for Yoshida's theorem hold, the weight of a weighted- 
homogeneous integral of (37) should be one of the resonances (42). Indeed, we find an 
integral by singularity analysis for any of the positive resonances given by equation (42). The 
integrals are given by the recurrent relations: 



m/2-1 



(2*0 

n— m— l+2k 



n 



(43) 



fc=0 



j=n— m+l+2k 



where the convention 1^ = 1 and 1^ = for all m, n with n = 2,3,... and m > n is 
adopted. 

For n = 3, the resonances are r = —3, —1,2 and a polynomial integral is 



r(2) 



C 2 = V>2 + ti| + 2uiU2 + 2?i2li3 



(44) 



However, it is simple to see that this is not the only first integral of the system (37). Defining 
u = U2 + U3, and using any constant value c of the integral (44) , the equations of motion take 
the form 

(45) 



2 V ; 



Integration of (45) yields 



u = — ccoth( 



c(t-to)- 



Using u instead of U2 as a new independent variable yields the equation: 

n 3 + uu 3 = U3 



(46) 
(47) 
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Which can be solved for u% yielding 



u , = _ c CQ S h(c(t-t ))-l 

sinh(c(i - t )) - c{t - to) + 2c7 y J 

where 7 is an integration constant. By eliminating the time between the solutions (46) and 
(48), a new first integral of the original equations is found: 

1 l*2 + « 3 -C 2 Ul +u 2 + u 3 

I tr = -2 7 = - ln( ■ — ) + (49) 

c u 2 + n 3 + c U1U3 

which is a transcendental function of the variables Uj. This integral could not have been 
found by the initial polynomial ansatz. 



4 Conclusions 

In this paper we have explored the question of whether it is possible to recover the first 
integrals of systems of first-order nonlinear ordinary different equations involving quasi- 
polynomial functions of the independent variables based on the information provided by 
singularity analysis. The main conclusions are: 

a) The theorem of Yoshida (1983) constrains the choice of an initial ansatz for an integral 
with undetermined parameters, leaving, however, an infinity of possible choices. 

b) The condition of Yoshida's theorem (V<3?(6j) 7^ and finite) holds for all the integrable 
Hamiltonian systems of two degrees of freedom included in Hietarinta's (1983) table, if hi is 
set equal to the principal balance 6j / for alH = 1, . . . 4. Other types of balances have to 
be considered in more general systems. 

c) Substitution of the Painleve series in a quasi-polynomial function $(x; c), where c is the 
vector of undetermined parameters, allows to separate the terms in powers of the time and 
determine the parameters c by singular value decomposition. Thus the information on quasi- 
polynomial integrals is contained in the Painleve-type series solutions around any movable 
singularity. 

d) In the case of balances r A with A > 0, left Painleve series must be used in the imple- 
mentation of the algorithm. 
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